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ABSTRACT 

I show that the temperature of nuclear star clusters, starburst clusters in M82, compact high- 
z galaxies, and some globular clusters of the Galaxy likely exceeded the ice line tempera- 
ture (Ti ce ps 150 — 170 K) during formation for a time comparable to the planet formation 
timescale. The protoplanetary disks within these systems will thus not have an ice line, de- 
creasing the total material available for building protoplanetary embryos, inhibiting the for- 
mation of gas- and ice-giants if they form by core accretion, and prohibiting habitability. 
Planet formation by gravitational instability is similarly suppressed because Toomre's Q > 1 
in all but the most massive disks. I discuss these results in the context of the observed lack 
of planets in 47 Tuc. I predict that a similar search for planets in the globular cluster NGC 
6366 ([Fc/H] = —0.82) should yield detections, whereas (counterintuitively) the relatively 
metal-rich globular clusters NGC 6440, 6441, and 6388 should be devoid of giant plan- 
ets. The characteristic stellar surface density above which T\ cc is exceeded in star clusters 

is ~ 6 x 10 3 M Q pc -2 r 1 '^, where /d gl mw is the dust-to-gas ratio of the embedding 
material, normalized to the Milky Way value. Simple estimates suggest that ~ 5 — 50% of 
the stars in the universe formed in an environment exceeding this surface density. Caveats and 
uncertainties are detailed. 

Key words: galaxies: formation — galaxies: evolution — galaxies: starburst — galaxies: star 
clusters: general — planets and satellites: formation — protoplanetary disks 



1 INTRODUCTION 

A key question is whether planet formation over cosmic time is 
concomitant with star formation, and whether or not the efficiency 
of planet formation is independent of environment and formation 
epoch. The correlation observed between giant planet frequency 
and host star metallicity by Gonzalez (19970 indicates that it is 
not, and implies that the z ~ planetary mass density Q p is at 
minimum a convolution of the enrichment history of the ISM of 
galaxies and their star formation rates. 

However, the relationship between fi p and the z — stellar 
mass density SI* is likely more complicated than a simple function 
of metallicity. Much of the star formation in normal spiral galax- 
ies, starbursts, and the rapidly star-forming galaxies at high red- 
shift occurs in massive star clusters (Lada & Lada 2003; Portegies- 
Zwart et al. 2010; Kruijssen 2012). The dense environment ac- 
companying stellar birth in these systems is likely to significantly 
impact the physics of planet formation. For example, dynamical 
interactions, UV photoevaporation, and superaovae have all been 



1 See also Laughlin 2000, Gonzalez et al. 2001, Laws et al. 2003, Santos 
et al. 2001, 2004, and Fischer & Valenti 2005. 



suggested as mechanisms that might inhibit, disrupt, or affect the 
formation, migration, and survival of planets in dense stellar en- 
vironments (Laughlin & Adams 1998; Armitage 2000; Adams et 
al. 2004; Fregeau et al. 2006; Adams et al. 2006; Fatuzzo & Adams 
2008; Gorti & Hollenbach 2009; Gorti et al. 2009; Proszkow & 
Adams 2009; Spurzem et al. 2009; Olczak et al. 2012). 

Because the timescale for these processes to affect planet 
formation is comparable to the planet formation timescale itself 
(~ Myr), and because most stars form in clusters, the planet forma- 
tion efficiency is likely a function of the stellar density of the birth 
environment. Since the bulk of the universe's stars are formed at 
z ~ 1 — 2, where the conditions for producing dense star clusters 
with high surface densities were common, but metallicities were 
generally lower, the functional dependence of Q p /Q ic (z) may be 
complicated. 

These issues are of relevance now because observations over 
the last decade have begun to constrain the incidence of planets in 
the local field disk population, in both open and globular clusters 
of the Galaxy, and, through microlensing, the incidence of planets 
in the Milky Way's bulge. These various data sets and future exten- 
sions have provided and will provide a wealth of information on the 
(likely) multi-dimensional space of planet formation efficiency as a 
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function of birth environment. In particular, Gilliland et al. (2000) 
and Weldrake et al. (2005) find a statistically significant lack of 
"hot Jupiters" in the globular cluster 47 Tuc. A multitude of tran- 
sit surveys in open clusters have begun to constrain the efficiency 
of hot Jupiter formation in these less dense systems (Mochejska et 
al. 2002, 2004, 2005, 2006; Hartman et al. 2005; Hidas et al. 2005; 
Rosvick & Robb 2006; Burke et al. 2006; Montalto et al. 2007; 
Hartman et al. 2009; Nascimbeni et al. 2012; see the analysis in Van 
Saders & Gaudi 201 1), and RV surveys have provided tighter con- 
straints on the population and, recently, some detections (Cochran 
et al. 2002; Paulson et al. 2004; Quinn et al. 2012; Pasquini et 
al. 2012). As I emphasize here and throughout this paper, such sur- 
veys, if carried out for a number of globular and open clusters with 
a variety of metallicities and stellar densities, might provide the first 
quantitative measure of the planet formation efficiency as a func- 
tion of both metallicity and stellar surface density. 

Complimentary to transit and RV surveys is microlensing, 
which could in principle measure the planet frequency in the Galac- 
tic bulge via self-lensing (Kriaga & Paczynski 1994) relative to the 
frequency in the foreground disk (e.g., Han & Gould 1995). Mi- 
crolensing may also be used to constrain the planet frequency via 
self-lensing in M31 (Baltz & Gondolo 2001; Chung et al. 2006). 
Since the Galactic bulge has had a different formation and enrich- 
ment history than the disk or halo, one would expect a different 
planet formation efficiency per unit star formed. Studies of the low 
mass IMF in the Galaxy's bulge and disk reveal a large difference 
in slope (see Fig. 8 from Zoccali et al. 2000) that might manifest 
an underlying difference in the relative planet population. Future 
observations of the relative planet frequency in these two envi- 
ronments (bulge and disk) might be interpretable with constraints 
on the planet population in the many different open and globular 
clusters of the Galaxy, again with a range of metallicities and sur- 
face densities. These studies would be further complemented by 
searches for the NIR signatures of protoplanetary disks in young 
and dense stellar environments, such as in Arches and NGC 3603 
(Stolte etal. 2004, 2010). 

With this backdrop, here I consider the importance of the 
strong infrared irradiation expected in a newly-born star cluster. 
Different from works focusing on photo-evaporation of protoplan- 
etary disks, I show that the UV and optical emission from massive 
stars, which is efficiently absorbed and scattered by dust grains, 
creates a thermal bath similar to a hohlraum. This strong thermal 
cluster irradiation is the focus of this work. 

I show that the heating and thermal structure of protoplanetary 
disks around Solar-type stars during the first few Myr of cluster 
evolution is in fact dominated by the irradiation from the cluster as 
a whole, and not from the host star, for semi-major axes larger than 
a few AU. This keeps the disks essentially isothermal and hot as 
a function of radius, affecting their structure and their potential to 
form planets via gravitational instability or core accretion. In fact, 
I show that in the interior of massive embedded star clusters, the 
temperature approaches and exceeds the condensation temperature 
of water ice, Ti cc ~ 150 - 170 K (Podolak & Zucker 2004; Lecar 
et al. 2006). Because the disks equilibrate to midplane temperatures 
at or above the value corresponding to the incident flux, this means 
that protoplanetary disks in such environments have no ice line. 
Stated succinctly, the central temperatures of many star clusters in 
formation exceed the fundamental temperature scale for planet for- 
mation by core-accretion: Ti cc . This decreases the total amount of 
condensable material by a factor of ~ 3 — 5, and thus inhibits giant 
planet formation (Lodders 2003; Lecar et al. 2006). Over a wider 



and less extreme range of cluster properties, the temperature is still 
very high, and this should in general suppress planet formation by 
gravitational instability. Even more generally, I find that in the first 
~ 1 — 10 Myr of evolution that the effective temperature of clus- 
ters and their embedded disks should be large than many tens of 
Kelvin, and for this reason studies of protoplanetary disk evolution 
for a much wider range of irradiation than has previously been con- 
sidered should be carried out. 

In Section[2] I estimate the temperature of embedded star clus- 
ters, and in Section [3] I estimate their protoplanetary disk temper- 
atures. I then compare the relative importance of cluster irradia- 
tion and (1) host star irradiation in passive disks, and (2) accre- 
tion power in active disks. Generically, I find that cluster irradiation 
dominates the thermodynamics and structure of passive and active 
disks outside a semi-major axis of ~ 1 — 5 AU around Solar-type 
stars. In Section[4] I discuss planet formation via core accretion and 
gravitational instabilty, focusing on the physics of the ice-line and 
the Q-criterion. I find that for star clusters above a stellar surface 
density of ~ 5 x 10 3 Mq pc -2 there is no ice-line and that, for 
even modest surface density clusters, the Q-criterion indicates that 
only protoplanetary disks with masses larger than ~ 0.2 Mq can 
be gravitationally unstable around Solar-type stars. 

In Section [5] I estimate the temperature at formation for a 
variety of systems, including Galactic star clusters and globular 
clusters, super star clusters in the nearby starburst M82, the cen- 
tral nuclear clusters of galaxies, and compact galaxies at high red- 
shift. Several globular clusters are likely to have exceeded Ti ce in 
at least part of their volume during formation, including 47 Tuc. I 
extend the discussion to continuously star-forming disks and dis- 
cuss planet formation in the context of local and high-z starbursts, 
rapidly star-forming galaxies, and AGN, finding that above a gas 
surface density of ~ 10 4 Mq pc -2 (or star formation surface den- 
sity of ~ 10 3 Mq yr _1 kpc~ 2 ) the temperatures of rapidly star- 
forming galaxies should exceed or approach Ti cc . Applying these 
results to the compact passive high-z galaxies, I find that most of 
the systems likely exceeded Ti cc during formation. 

In Section [6] I discuss the results, focusing first on uncer- 
tainties in the estimates, including variations in the IMF, radiation 
transport with distributed sources, estimates of the optical depth in 
star clusters, and the time dependence of cluster and planet for- 
mation. With these uncertainties in mind, I discuss local globu- 
lar clusters, and in particular suggest other transit/RV searches for 
hot Jupiters in Galactic globular clusters than have central temper- 
atures below Tico, but metallicities comparable to 47 Tuc. I pro- 
vide a sketch of a calculation of the fraction of all star formation in 
the universe that occurred above a given stellar surface density and 
metallicity, with implications for the z — planetary mass density. 
Section|7]provides a conclusion. 



2 THE TEMPERATURE OF STAR CLUSTERS 

A star cluster has mass M, radius R, and a stellar surface density 
£ = M /47r7? 2 If the stars form in a time less than 3 MyrQ then 



2 We define E with the "4" since we are interested in the total flux, as seen 
from within the cluster. 

3 The dynamical timescale is ~ 10 5 yr (.R/PC/E4) 1 ' 2 . 
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Figure 1. Effective temperature of star clusters as a function of time for an 
instantaneous burst of star formation for E = 10 5 — 10 2 Mq pc — 2 (top 
to bottom). The gray band shows Ti cc ~ 150 - 170 K. After t ~ 3 Myr, 
with the normalization set by eq. (|2j- 



the flux at formation can be estimated as 



F ss *E ~ 1.6 x 10 7 L 



PC 



vi/ 



1500 ^4, 



(1) 



where E x = E/IO^Mq pc~ 2 and Vfisoo = (tf/1500) L /M 
is the light-to-mass ratio for a zero-age main sequence (ZAMS) 
stellar population assuming a fully-populated standard IMF (e.g., 
Leitherer et al. 1999)Q This flux corresponds to an effective tem- 
perature of 

vl/4 



T eff ~ 101 K (*1500 E 4 ) 



(2) 



Figure[T]shows T e s as a function of time for an instantaneous burst 
of star formation, using a Salpeter IMF from 1 — 100 Mq, for 
several values of E. After 3 — 4 Myr, and over the first 100 Myr 
the time-dependence of T e g(t) is well approximated by f~ 1//3 . 
The gray band shows the temperature for water ice sublimation, 
lice =± 150 - 170 K. The effect of variation in the IMF on * and 
Teff is discussed in Section RTTl 



Since the focus here is on the application to protoplanetary 
disks, one can ask about the total flux absorbed by each side of 
an optically-thick disk situated at the center of an optically-thin 
spherically-symmetric star cluster with radial emissivity profile 
j(r) — Vl/p. For illustration, taking a Hubble profile for the stel- 
lar density distribution, p(r) = po/[l + (r/ro) 2 ] 3 / 2 , one obtains 
F(r = 0) = *po^o/4 = ctsbTcV implying 



Teff ~ 108 K (*1500 po,5 ni.pc) 



1/4 



(3) 



where ro iPC = ro/pc and po,5 = Po/1O 5 M0 pc 3 so that 
the equivalent stellar surface density inside ro is E = M(r < 



ro)/47rro — 1.7 x 10 4 Mq pc -2 . More generally, the total energy 
density at any radius r within an optically-thin cluster of emissivity 
j(r) = S&p(r) is 



u(r) = 1 



j(r') L- In 



1 + r'/r 



1 — r'/r 



dr' . 



(4) 



2.1 Optically-Thick Clusters 

If the star cluster is compact compared to the optically-thick em- 
bedding dust/gas distribution, the temperature of the star cluster is 
larger than T c ff because the radiation field is trapped and diffusive, 
as in a simple stellar atmosphere. In this case, the cluster tempera- 
ture is 



r 4 

-*■ r. 



TR + 



T 4 



(5) 



where tr is the Rosseland-mean optical depth through the sur- 
rounding dusty medium. For the purposes of analytic estimates, in 
the temperature regime 160 K< T < 1500 K, I approximate the 
Rosseland-mean dust opacity as constant: 



5/dg, mw cm 2 g 1 , 



(6) 



where I assume a Galactic dust-to-gas ratio (/dg, mw = !)■ Com- 
paring with Adams & Shu (1986), Bell & Lin (1994), and Semenov 
et al. (2003), this approximation agrees with more detailed calcu- 
lations to a factor of ~ 2 over this temperature range. For lower 
temperatures (T < 160 K), 



KR ~ Ko/dg,MwT , 



(7) 



where k ~ 2.4 x 10" 4 cm 2 g _1 KT 2 (Bell & Lin 1994; Semenov 
et al. 2003). 

Approximating the cluster as a geometrically thin shell of gas 
surface density E 9 and gas fraction f 9 = E 9 /E, the total optical 
depth is 



TR ~ Kfl/gE ~2Kfl /9E4. 



(8) 



For T c > 160 K and tr ^> 2/3, the cluster temperature is then 
approximately 



T c ~ 170 K (/ dg , mw U *i5oo) 1/4 E 1 / 2 , 
while for lower temperatures (T < 160 K)[f| 

T c ~ 130K(/ dg ,Mw/ 9 *i50o) 1/2 E 3 .8. 



(9) 



(10) 



The radiation energy density is u = a r T 4 , where a r is the radiation 
constant. These estimates for hot, high surface density embedded 
clusters are similar to those derived by Murray (2009). 



2.2 Optically-Thin Clusters 

When the cluster is optically-thick to the UV/optical radiation from 
the stars, but optically-thin to the re-radiated FIR emission, the tem- 
perature is T 4 ~ T* b /Atr, and the energy density is given by 
equation ©■ The critical surface density below which the cluster 



4 For clusters with small total mass < 10 3 Mq, the assumption of a fully 
populated IMF breaks down and \t is decreased. 



5 There is only a small dynamic range in f g S where tr > 2/3 and T c < 
160 K. See Figured 
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Figure 2. Temperature as a function of cluster stellar surface density. The 
thick solid black line is T c g for the cluster (eq. |2) for t < 3 — 4 Myr 
(see Fig. [TJ. The red lines show the central cluster temperature T c for 
/dg, MW x fg = 100, 10, 1, 0.1, and 0.01 (eq.[5j, using the cluster optical 
depth, tr, and using the approximate opacities of eqs.|6]and|2] A numerical 
solution to eq. using the full opacity tables of Semenov et al. (2003) for 
/dg, MW = / 9 = 1 is shown by the dotted line. A protoplanetary disk 
that is optically-thick to the cluster's FIR emission at temperature T c at- 
tains a minimum temperature T mid ~ T c if the cluster has tr > 1 . For 
tr < 1, a protoplanetary disk that is optically-thick to the cluster's FIR 
emission at temperature T c attains a minimum temperature T m ;d ~ T e ff . 
The gray band shows Ti cc ks 150 — 170 K. The horizontal dashed lines 
show condensation temperatures for troilite, magnetite, ammonia, methane, 
and argon ices (highest to lowest; Lodders 2003). 

has tr < 2/3 can be derived roughly using equations ((5} and ((7} 
as E TR = 2/3 ~ 1500 M Q pc" 2 (/ dg , «w/ 9 )- 2/3 C. 



2.3 UV & Optically Optically-Thin Clusters 

The cluster is optically-thin to both the UV/optical radiation from 
massive stars and to the re-radiated FIR emission when E < 

5 Mq pc -2 (/dg, mw/s ) _ 1 > which may occur after the cluster has 
evolved out of the embedded phase. In this limit, T e g is still given 
by equation (O and u(r) is given by equation but the tempera- 
ture of the radiation field is of course much hotter (~ 10 4 K) than 
in the limits of Sections |2. 1 1 and [2~2l This hard radiation field will 
super-heat dust grains at the surface of protoplanetary disks to a 
temperature larger than T e ff because of the emissivity properties of 
grains in the UV and FIR. 



3 THE TEMPERATURE OF PROTOPLANETARY DISKS 

Given the ambient radiation field in the cluster, one can estimate its 
importance for the thermal structure of protoplanetary disks around 



individual stars within the cluster, and compare the importance of 
cluster irradiation to the insolation from the host star and to the 
thermal structure determined by accretion in active disks. 

3.1 Embedded Disks 

There are three primary regimes of interest. To illustrate these, I 
adopt a simple model of a protoplanetary disk around a star of mass 
M^, radius 71*, temperature T*. The disk has gas surface density 
E = Eo(AU/a) a: lf| where a is the semi-major axis, disk midplane 
temperature T m id set by radiative and hydrostatic equilibrium, and 
vertical optical depth Ty — K_fj(r m id)E/2, where kr is given by 
equations © and {7} if the grain size distribution and composition 
within the disk is the same as within the star cluster and if we con- 
sider only temperatures below the dust sublimation temperature at 
~ 1500 K0 

Neglecting all other sources of disk heating, including that 
from the central star and accretion, the three regimes of clus- 
ter irradiation for the disk are as follows. At small a the disk 
is optically-thick to radiation of temperature T c , and to its own 
radiation T m i^. Here, T m id ~ (u/a r ) 1//4 . Next, at an interme- 
diate range of a, the disk may be optically-thick to radiation of 
temperature T c , but optically-thin to its own radiation T m id and 
Tmid ~ (u/a r Tv) 1//4 . Finally, at large a, the disk is optically-thin 
to radiation of temperature T c and to its own radiation T m id and 
T mid ~ [e(T c )/e(T mid )] 1/4 (u/a r ) 1/4 , where e is the dust emis- 
sivity at temperature T[f| 

These limits show that a protoplanetary disk embedded in a 
cluster of effective temperature T e g (eq. [2]l has ~ T e g as its min- 
imum temperature. If the cluster is optically-thick, then the mini- 
mum temperature of the optically-thick regions of the disk is ~ T c 
(eqs.|9l&|10ll, which exceeds T e g by a factor of ~ r 1 ^ 4 (eq.[5). 

Figure [2] shows several temperatures of interest as a func- 
tion of cluster stellar surface density for a newly-born cluster 
(t < 3 - 4 Myr). The thick solid black line shows T cB (eq. O 
for * = 1500L Q /M Q (see Fig.QJ. The red solid lines show T c for 
tr > 2/3 obtained by solving equation l[5), using the simplified 
opacities of equations (f6]l and Q, and for /d g , mw X f g = 100, 10, 
1, 0.1, and 0.01. The break in each red line occurs at T c — 160 K at 
the break in the assumed opacity law (eqs.|9l&|10t. A more com- 
plete treatment of the Rosseland-mean opacity, as presented in Bell 
& Lin (1994) or Semenov et al. (2003), introduces inflections in 
T C (E) at each break in the opacity curve, as shown in the dotted 
curve, which gives the full solution to equation ((5} using the opac- 
ity tables of Semenov et al. (2003) for /dg, mw X f g = 1. The 
close correspondence between the dotted curve and the red curve 
validates our use of equations ((6) and ((7} in making analytic esti- 
mates below. The gray band shows the temperature for water ice 

6 For a minimum mass Solar nebula, Eo ~ 10 3 g cm -2 and x ~ 3/2 
(Weidenschilling 1977). 

7 Note that the grain size distribution of proto-planetary disks is likely dif- 
ferent from that of the ISM, with a larger maximum grain size and different 
slope for the power-law size distribution (e.g., D'Alessio et al. 1999, 2001). 
These differences quantitatively affect the dividing lines between the differ- 
ent semi-major axis regimes discussed below, but all qualitative conclusions 
remain unchanged. More discussion of this issue is provided in Sections [3.3| 
andH4l 

8 Note that t(T c )/e (T m ;d ) would only be much different than unity in the 
case of a UV/optically optically-thin cluster as discussed in Section l231 
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sublimation Ti ce — 150 — 170 K. The short horizontal dashed lines 
show condensation temperatures for troilite (~ 700 K) and mag- 
netite (~ 370 K), and then the ammonia, methane, and argon ices 
(Lodders 2003). For times > 3 — 4Myr, T e g and T c decrease ap- 
proximately as i~ 1//3 (Fig.[T|l. 



3.2 Cluster Irradiation versus Host Star Irradiation 

It is useful to compare the incoming flux from the combined stellar 
population of the star cluster with the flux from the central star in a 
passive flared disk (Chiang & Goldreich 1997; hereafter CG97): 



4a 2 



1/4 



-3/7 
'10 AU 



4 /2 ' 



(id 



where a ~ 0.005(17^ + 0.05a 2 /^ is the flaring angle, Tq = 
T*/6000K, Rq = R^/Rq, and where the second approximate 
equality follows from ignoring the first term in a at large a. 

As discussed in CG97, the hot stellar irradiation produces a 
superheated dust layer (of temperature Tdust) at the surface of the 
disk, which re-radiates approximately half of the absorbed light 
back into the disk, maintaining T m id. As in the case of cluster ir- 
radiation alone, there are three regimes of irradiation from the cen- 
tral star separated by semi-major axis: (1) at small a (< 100 AU 
for CG97 parameters), the disk is optically-thick to both Td us t and 
Tmid; (2) at an intermediate range of a the disk is optically-thick 
to Td ua t, but optically- thin to radiation of temperature T m id; (3) at 
large a (> 200 AU in CG97) the disk is optically-thin to both. In 
all cases, T m jd for the passive disk is set by the incoming flux from 
the central star, T e fi, * in equation dl II ). but with an overall correc- 
tion for the fact that only half of the incoming flux is radiated into 
the disk, a correction for the vertical disk optical depth in regime 
(2), and a correction for the ratio of emissivities between Tdust and 
T mid in regime (3). 

We can now compare the importance of cluster irradiation to 
that from the central star directly. When the cluster is optically- 
thick to its own FIR radiation (tr > 2/3), the relevant comparison 
is between T c in equations ((9} and d 10b and T c s, *. Setting T c = 
T c ff , *, one derives the critical semi-major axis outside of which the 
cluster radiation field dominates the radiation field from the central 
star: 

7/12 

(T c > 160 K) (12) 



ttcrit 



0.9 AU 



Rt,T r 



dg, MW 



fg *1500S| 



0.6 AU 



RqTq 



/dg, MW fg #1500 £4 



7/6 



(T C <160K) (13) 



When the cluster is optically-thin to its own FIR radiation the rele- 
vant comparison is between T g from the cluster (eq.|2} and T e ff : *: 



3.0 AU 



0'0 



4 n7/12 



#15 



{tr < 2/3). 



(14) 



Since at these small values of the semi-major axis the disk is 
optically-thick to the incoming FIR radiation from the cluster, the 
re-radiated emission from the hot surface disk dust layer, and to 
the re-radiated emission from the disk midplane, these comparisons 
show that the cluster radiation field dominates the temperature and 
structure of passive irradiated disks for a > 1 AU around typical 

stars in a cluster with E = 10 4 Mo pc -2 ; for a < a cr it, T m id is 



1000 



500 



too 




100 



Figure 3. Protoplanetary disk temperature as a function of semi-major axis 
a for a star with M* =Mq, T* = T = 6000 K, R+ =Rq, and 
a disk with £ = 3 X 10 3 g cm -2 (AU/a) 3 / 2 . The gray band shows 
Tico = 150 - 170 K, the short dashed black line shows Tcff^/2 1 / 4 
(eg. I UK which approximates the contribution to the midplane temperature 
from stellar irradiation alone. The dotted line shows the analytic estimates 
°f T m id, acc assuming an accretion rate of 10 -8 Mq yr _1 using eqs. 1 1 6t 
and i ni . The heavy black lines show the calculated value of midplane tem- 
perature of the disk as a function of a, taking into account host star in- 
solation, accretion, and cluster irradiation using the opacities Semenov et 
al. (2003) with T c = 50, 100, and 200 K, respectively, from lowest to 
highest. Higher T c in general moves the ice-line to larger a, and when 
T c > Ti ce (e.g., forT c = 200 K), the ice-line disappears entirely. 



set by stellar irradiation, but for a > et C rit , the cluster sets T m id. Its 
minimum value is T c g (black solid line in Fig.|2j in the case of an 
optically-thin cluster, and it reaches T c in the case of a cluster that 
is optically-thick to its own FIR radiation. 

It is worth noting that because the geometry of cluster irradi- 
ation is different than for the central star, the vertical structure of 
the disk will be modified. In particular, the flaring angle in a pas- 
sive disk irradiated only by a central star is expected to scale as 
h/a oc a 2//7 (CG97), whereas for a passive isothermal disk irra- 
diated isotropically by the cluster one expects h/a oc a 1 / 2 . This 
stronger flaring will increase the insolation from the host star, and 
in a self-consistent model for the disk, cluster, and star, and one 
expects the estimates above for a cr i t to be modified quantitatively. 



3.3 Cluster Irradiation versus Accretion 

For an active disk with steady accretion rate M-g = M /10 -8 Mq 
yr _1 around a star of mass M* measured in units of Mq, the flux 
from the disk is 



T 



off, acc 



3 MQ, 2 



-i 1/4 



87T <7SB 



26K (M_ 8 M t /a 3 6AU ) 1/4 . (15) 
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Ignoring other sources of heating, the midplane temperature of the 
disk is larger than T e s l acc by a factor of [(3/4)rv] 1 ^ 4 . Taking val- 
ues characteristic of a minimum mass Solar nebula for illustration 
— E(a) = E (AU/a) 3/2 with E ,3 = E /10 3 g cm" 2 — the 
midplane temperature is 



Tmid.acc « 260K(E , 3 M_ 8 M,) 1/4 a 



9/8 
'2 AU 



(T mid > 160 K) 
(16) 



T 



mid, acc 



59K(E , 3 M_ 8 M*) 



1/2-9/4 



a 5 AU (^mid < 160 K), 



(17) 

where the top expression assumes equation ((6} and the lower ex- 
pression assumes equation 00 The above scalings show that ac- 
cretion heating generally dominates T m id at a < few AU, whereas 
passive heating from the central star dominates T m id at larger semi- 
major axis (CG97; compare eqs.|17|andll It. 

In order to evaluate the importance of cluster irradiation for 
the structure of active disks, the relevant comparison is then be- 
tween T m id, acc and either T c in the case of a FIR optically-thick 
cluster or T e g in the case of a FIR optically-thin cluster. Roughly, 
T m id wll l on ly be strongly modified when max[l, trJ'I'E for the 
cluster exceeds TyMfl /8 for the optically-thick regions of the 
protoplanetary disk. For optically-thick clusters with tr > 2/3 the 
critical semi-major axis beyond which cluster irradiation dominates 
accretion heating in setting T m id is 



2.9 AU 



E , 3 M_ 8 M* 



-, 2/9 



(18) 



_/dg, MW fg Vf r 150oE 4 

which is correct for either T c , T m id > 160 K or T c , T m id < 
160 K. When the cluster is optically-thin (tr < 2/3), equations 
© and {FT} imply that 



fflcrit 



3.9 AU 



(E , 3 M_ 8 M*) 2 



^15(X)E4 



1/9 



(tr < 2/3). 



(19) 



These comparisons show that the cluster radiation field dom- 
inates the disk structure of active accreting disks for semi-major 
axes larger than ~ few — 10 AU around a typical star with an ac- 
tive accretion disk embedded in a cluster with E > 10 2 Mq pc -2 . 

The weak dependencies of a cr it on E and M in equation ( 1 1 9b guar- 
antee that beyond a ~ few — 10 AU, cluster irradiation dominates 
accretion, even for accretion rates as high as 5 x 1O _7 M0 yr~ 4 
and cluster surface density E as low as 10 2 ' 5 Mq pc -2 . 



However, as I discuss in Sections 14.21 and 14.31 the cluster- 
irradiated massive disks most prone to gravitational instability will 
have accretion rates as much as ~ 10 2 — 10 4 times higher than used 
in the estimates above, and particularly at early times. The normal- 
ization of equations j 1 8| i and d!9t increases to ~ 22 and ~ 11 AU, 
respectively, for M — 10~ Mq yr _1 . Additionally, the accretion 
luminosity in such a system would dominate the passive heating 
contribution in Section |T2l pushing the critical semi-major axis at 

9 For a proto-planetary disk dust grain size distribution, the Rosseland- 
mean opacity is decreased by a factor of ~ 2 — 3 relative to the estimates 
of equations |6} and j7} in the temperature range of interest (D'Alessio et 
al. 1999, 2001). In general, this lowers the midplane temperature of active 
disks with respect to the estimates in equations | |161 and (TTJ- This moves 
the ice-line in active disks to smaller semi-major axes (Sasselov & Lecar 
2000; Lecar et al. 2006). Setting equation (Tg) equal to T Icc = 170 K 
one finds that a Icc ~ 2.8 AU[(k r /4 cm 2 /g)E , 3 M_ 8 M] 8 / 9 , whereas 
a calculation with the D'Alessio et al. (2001) opacities by Lecar et al. (2006) 
finds ai cc ~ 1.7 AU for the same parameters. 



which the cluster irradiation dominates the host star to much larger 
values than given in equations d 1 21>-(l 1 4I> . Most importantly, such 
high accretion rates would significantly decrease the disk lifetime. 



3.4 A More Complete Calculation 

In order to make these comparisons between cluster irradiation and 
host star irradiation and accretion explicit, in Figure |3]l show the 
midplane temperature of a protoplanetary disk (thick solid lines), 
including all three effects, for T c = 200, 100, and 50 K (highest to 
lowest). 

The dashed line ("stellar irradiation") shows T c ff . * (eq. [TTJ 
CG97), and the dotted line ("accretion") shows the analytic ap- 
proximations for Tmid, acc in equations J 1 6b and | |17I >, The midplane 
temperature is calculated by solving the implicit equation (see, e.g., 
Sirko & Goodman 2003) 



4 

T V + - 



2 

3iy 



Tcff, acc + f 1 H ) (T 4 ff . * + T 4 ) , 

(20) 

where Ty(T m id) is the vertical optical depth to the disk's own ra- 
diation, and where I have assumed the opacities of Semenov et 
al. (2003). As in CG97, the last term allows for the fact that the disk 
may be optically-thick to the incoming radiation field, but optically- 
thin to its own re-radiated emission. 

The calculation is not self-consistent in the sense that it as- 
sumes M is constant throughout the disk, and that T c g, * is given by 
equation i ll It . In reality, the disk structure, and hence the amount of 
stellar irradiation as a function of a, should deviate from equation 
i ll It because of accretion, and because the disk becomes isother- 
mal (Tnid — ► T c ) at large a. In addition, because the scale height 
is determined without reference to M, in such a disk as plotted in 
Figure [3] the viscosity parameter a (not to be confused with a in 
eq.ll 111 would continuously vary. For these reasons, the calculations 
shown in Figure |3]are meant to be a sketch. Even so, Figure [^dra- 
matically illustrates the importance of cluster irradiation. Even for 
T c = 50 K, the thermal structure of the disk is completely modified 
from the expectation of passive stellar irradiation for semi-major 
axes larger than ~ 10 AU. Comparing T c = 50 K to T c = 100 K, 
one sees that the ice line moves out slightly in a and that the disk 
structure is modified from the expectation of the simply active disk 
beyond a > 4 AU. For T c = 200 K, the ice line completely dis- 
appears, and the structure of the disk changes noticeably beyond 
> 2 AU. 

Detailed self-consistent calculations of active disks with both 
cluster and host star irradiation are left to a future effort. 



3.5 Other Regimes 

Many more regimes than those listed in equations dl21 -dl4t and 
dl8b-l|19l>, and represented in Figure [3] could be enumerated. For 
example, when the cluster is UV/optically optically-thin and the 
cluster has the hard radiation field characteristic of a young stellar 
population (~ 10 4 K, Section [2~3t an active and irradiated proto- 
planetary disk will be optically-thick to this incoming UV radia- 
tion at much larger semi-major axis than the FIR radiation of an 
optically-thick cluster, thus effecting at what semi-major axis the 
cluster, the accretion, or the insolation from the host star dominates 
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the disk structure and thermodynamics. In addition, such a clus- 
ter radiation field would produce super-heated dust layers on all its 
protoplanetary disks, thus affecting the dust-reprocessed radiation 
field the disks would see at larger a. 



4 PLANET FORMATION 

The estimates above show that the temperature within compact em- 
bedded star clusters is expected to be very high. Protoplanetary 
disks within these clusters can be expected to have a minimum tem- 
perature of order T e ff given by equation lO, and with time evolution 
shown in Figure [T] Highly embedded clusters reach considerably 
higher central temperatures, as shown by the red lines in Figure 
[2] for different values of the dust-to-gas ratio and cluster gas frac- 
tion. As discussed in Section[3]and shown in Figure[3] the radiation 
field from a young cluster is likely to dominate the thermodynam- 
ics and structure of both passive and accreting protoplanetary disks 
on scales larger than ~ 1 — 10 AU, depending on the parameters 
of the system in consideration. Here, I consider the zeroth-order 
effects on planet formation in the core accretion and gravitational 
instability pictures for the formation of gas- and ice-giant planets. 



4.1 Core-Accretion 

The dominant picture for the formation of gas and ice giant plan- 
ets is the core accretion model (Mizuno 1980; Pollack et al. 1996), 
which requires a protoplanetary core to grow in the disk to a criti- 
cal mass of ~ 10 M®, thereby initiating runaway accretion of gas 
from the disk and rapid growth to the isolation mass (e.g., Armitage 
2010). Critical to the core accretion model is the "ice line," the 
semi-major axis where the temperature of the disk drops below the 
sublimation temperature of water ice, because here the surface den- 
sity of all condensables is ~ 3 — 5 times higher than that in only 
refractories. It is this "extra" material, which exists only outside 
the ice line, that enables the formation of ~ 10 M® mass cores and 
the core accretion phenomenon. In standard models of disk evolu- 
tion including irradiation from the central star and accretion, the 
ice line typically occurs at a ~ 2 — 3 AU (e.g., Sasselov & Lecar 
2000; Lecar et al. 2006; Kennedy & Kenyon 2008), and it is this 
fundamental scale that is posited to explain the dichotomy between 
the terrestrial and gas/ice giants in the solar system. 

Looking at Figure [2] one finds that in star clusters with high 
surface densities and/or high gas fractions or dust-to-gas ratios, T c 
is greater than the ice line temperature of Ti ce — 150 — 170 K for a 
fully populated initial mass function, for times less than ~ 3 Myr. 
Setting T c = Ti cc in equation ((9}, one derives the critical stellar 
surface density above which T c > Ti ce and there is no ice line in 
the star cluster is 

E Icc ~ 6 x 10 3 M pc" 2 (Jj^) (/d g ,Mw/ 9 *i50o)" 1/2 . 

(21) 

Note the metallicity dependence in the denominator of this expres- 
sion, which shows that a high value of the dust-to-gas ratio drives 
the system to higher temperatures, which could completely inhibit 
giant planet formation in the core accretion theory. As shown in 
Figure [2] and equation &2l\ , for super-solar metallicities or larger 
gas fractions, Ei co decreases: for /a gl mw X f B — 100, Ei cc ~ 
4 x 10 2 M© pc -2 . In contrast, for much lower /d g , mw, as might 




1000 1 1 — 1 — 1 — 1 J 

100 1000 10" 10 5 

£ (M pc-») 

Figure 4. Critical protoplanetary surface density normalization Eo, cr it (see 
eqs. 1221241 above which the disk has Q < 1 for f 3 X /jg mw = 100 — 
0.01 (top to bottom; same as Fig. [2), assuming £(a) = S (AU/a) 3 / 2 . 
The normalization for the MMSN is shown for comparison as the gray band. 
Very large Sq , and correspondingly large disk mass, is required for Q < 1 
in stellar clusters of essentially any surface density. 

occur in metal-poor proto-globular clusters with [Fc/H] = —2 so 
that /d g , mw X f g = 0.01, Eico increases to ~ 5 x 10 4 Mq pc~ 2 . 

If giant planets in fact form by core accretion, then because 
this is a critical phenomenon requiring a large amount of condens- 
able material, one expects clusters bora with E > Ei co to have 
only terrestrial rocky planets, and for those planets to be devoid of 
water. Such a disk is sketched in Figure [3] which shows that for 
T c — 200 K (top thick solid line), the ice line simply disappears. 
I discuss the search for giant planets in the relatively metal-rich 
globular cluster 47 Tuc ([Fe/H] = —0.76), the discovery of a giant 
planet in NGC 6121, and other globular clusters in Sections [5] and 

m 



4.2 Gravitational Instability 

Gas- and ice-giant planets may also form by gravitational insta- 
bility (Boss 1997, 2000; Pickett et al. 2003; Durisen et al. 2007). 
Planet formation by gravitational instability puts two necessary 
physical requirements on the structure of the protoplanetary disk 
that must be satisfied simultaneously: (1) Toomre's Q must be less 
than ~ 1 and (2) the cooling timescale must less than the local or- 
bital timescale (Gammie 2001). Recent work has shown that these 
requirements can only be met at large semi-major axes around stars 
with very massive disks (Matzner & Levin 2005; Rafikov 2005). 
Such massive protoplanetary disks that become gravitationally un- 
stable at large a have also recently been argued to produce low- 
mass binary stellar companions instead of planets (Boley 2009; 
Kratter et al. 2010; see also Bonnell & Bate 1994, Bate 2000). 

Effects on Q by irradiation have been considered in the con- 
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text of IR re-radiation by protostellar envelopes by Matzner & 
Levin (2005), and in the context of disk models by D'Alessio et 
al. (1997), Cai et al. (2008), Rafikov (2009), Vorobyov & Basu 
(2010), Rice et al. (201 1), Kratter & Murray-Clay (201 1), and Zhu 
etal. (2012). 

The most obvious consequence of strong irradiation in the 
cluster environments investigated here is that Toomre's Q is larger 
at fixed protoplanetary disk surface density relative to a non- 
irradiated protoplanetary disk of the same properties. Thus, larger 
surface density is required to drive the disk gravitationally unsta- 
ble. Taking c s = ^/k B T/fi with /i = 2.3m p , tt = (GM/a 3 ) 1/2 , 

= M*/Mq, and assuming that E s (a) = E (a /a) 3/2 , the 
critical temperature above which Q = 1 = c s £l/irGY, g scales as 
Tq = i oc Eq/M* with no a-dependence. In this case, the critical 
normalization Eo, cr i t above which the disk will have Q < 1 is 

Eo.crit ^ 5.7 x 10 4 g cm" 2 M 1/2 (*i5ooE 4 ) 1/8 (22) 

for an optically-thin cluster, and it increases to 

Eo.crit ^ 7.8 x 10 4 g cm" 2 M 1/2 (/ dg , M w/ 9 *i50oS!) 1/s 

(23) 

and 

E 0iCrit ^ 1.1 x 10 s g cm" 2 M 1/2 (/ dg ,Mw/ 9 *i50oSl) 1/4 

(24) 

for Tc > and T c < 160 K, respectively. These values are ~ 50 — 
100 times the typical value assumed for the MMSN (So = 1.5 — 
3 x 10 3 g cm -2 ). The implied disk masses are M D ~ 0.6-0.8 M Q 
for an outer disk semi-major axis of 50 AU. 

Figure [4] shows Eo, crit as a function of the cluster stellar 
surface density E for the same parameters as in Figure [2] f g x 
/d g ,MW = 100, 10, 1, 0.1, 0.01 (red lines, top to bottom; eqs.l23l 
124b and for T — T er r (black line; eq. 122b . The normalization for 
the MMSN is shown as the gray band. Thus, for disks embedded 
in star clusters with surface density E, and given f g x /d g , mw, 
the protoplanetary disks must have Eo > Eo, crit for Q < 1. This 
simple comparison assumes that T m id = T c , which is justified for 
large values of Eo, crit required by equations (122b- d24b since the 
vertical optical depth in the disk is larger than unity all the way out 
to ~ 100 AU. 

For different protoplanetary disk profiles, the scalings become 
more complicated. For E 9 (a) = T,o(ao/a) x , x = 1, and ao = 
AU, the critical temperature above which Q = 1 is 

T Q= i = 62 Ka 2 oAu Sg, 4 /M*. (25) 

I have purposefully scaled all quantities for a massive disk, with 
total mass Md ~ 27ra Eo x 20AU ~ 0.14E , 4M . Comparing 
with Figures[T]and|2]one sees that for most of the parameter regime 
considered T c g for the cluster exceeds Tq=i and that the disks will 
have Q > 1 for a < 20 AU. For the optically-thick clusters with 
T c > Tcff, T c can in many cases dramatically exceed Tq=i, as 
shown in Figure|2] 

However, note that Tq = i depends strongly on Eo, a, and M*. 
For an isothermal disk with x = 1, Q oc a _1//2 , and thus equation 
( 125 b can be interpreted as saying that in order for the disk to have 
a region with Q < 1 in a cluster with temperature equal to Tq-i, 
the disk need only extend to a semi-major axis beyond 20 AU (for 
the parameters of eq. 125 b . For example, equating Tq=i with the 
expression with the effective temperature of a star cluster of surface 
density E in equation {2^, one derives the critical semi-major axis 
of the protoplanetary disk beyond which Q < 1 in a disk with 



x = 1: 

a Q=1 ^33 AUA/ ^ 2 150oE 2 4)1/4 . (26) 

^0,4 a 0, AU 

Alternatively, one can think of this as a limit on the disk mass Md'- 

M d ,q = i ~ 0.18 M (M* a 20 Au) 1/2 (*i 5 ooE 4 ) 1/8 , (27) 

which says that for Md > Md, q=i in a cluster of stellar surface 
density E, the disk has Q < 1 at a ^ 20 AU. Thus, for gravitational 
instability to operate at, say, a = 100 AU in a disk with x — 1 
around a solar-mass star, in a cluster of nearly any stellar surface 
density E (since Md oc E 1 ' 8 ), the disk mass must be very large: 
M D ^ 0.2 - O.4M . 

The above estimates are only for optically-thin clusters. For 
the case of an optically-thick cluster with tr > 2/3, the limit on 
Md increases, but again depends on x. 

Overall, these estimates imply that cluster irradiation strongly 
suppresses gravitational instability by increasing Q above its nom- 
inal value for a disk in isolation. The additional requirement on the 
formation of planets by gravitational instability is that the cooling 
time of the disk should be less than the characteristic orbital time. 
The reader is referred to Kratter et al. (2009), Cai et al. (2008), 
Rafikov (2009), Rice et al. (2011), Kratter & Murray-Clay (2011), 
and Zhu et al. (2012) for a discussion of this criterion including 
irradiation. I note here only that the magnitude of the cluster irradi- 
ation indicated by Figures (QJ, and ((2), and which I estimate for ob- 
served stellar systems in Section|5] in general dramatically exceeds 
what has been so far considered in the literature for protoplanetary 
systems forming by gravitational instability. 



4.3 Accretion & Migration 

Because of the increase in sound speed at large semi-major axis in 
disks with strong cluster irradiation with respect to those without, 
one generically expects active disks to accrete at higher rates, yield- 
ing a shorter overall planet formation timescale, and for planets em- 
bedded within them to migrate faster, if the predominant migration 
mechanism is Type II. This conclusion follows from the accretion 
timescale in a disk with viscosity coefficient a: 

tvisc ~ ~ 4 x 10 4 yr c*o 1 (M* aiooAu) 1/2 TiooK- (28) 

a c| 

One also expects a somewhat larger accreting (active) atmosphere 
in layered disks with a dead zone (Gammie 1996). Similarly, in a 
self-gravitating marginally stable disk, the characteristic accretion 
rate is 

o 3 

M ~ ~ 10~ 4 M yr- 1 !^, (29) 

which also contains a strong temperature dependence. 

In addition, for disks embedded in clusters that exceed the ice 
line temperature one expects less total mass in planetesimals, and 
one then expects planetesimal migration to be less effective overall. 
This implies that the rocky planets that form in clusters exceeding 
the ice line temperature might have different migration history than 
their non-cluster analogs for at least three reasons: (1) there might 
be no Jupiter-mass planets to instigate, or regulate, planetesimal 
migration because the core-accretion process is quenched; (2) there 
will be less mass in planetesimals to begin with because ice cannot 
condense; (3) the viscous time for the disk will be short because 
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tvisc oc T" 1 , which is then proportional to the stellar surface den- 
sity of the cluster as t v i sc oc E~ 1//4 for an optically-thin cluster to 
tvisc oc E~ 1//2 for an optically-thick cluster with T c > Ti cc . These 
countervailing effects (1 + 2 versus 3) need to be considered in 
assessing the demographics of planets found (NGC 6121; see Sec- 
tions [6} or not found (47 Tuc) in dense stellar environments, 
even at fixed metallicity. 



5 COMPARISON WITH DATA & SPECIFIC SYSTEMS 

In this section I use the analytic estimates of the previous sections 
to estimate the temperatures of star clusters, globular clusters, and 
whole galaxies in formation. 



5.1 Super Star Clusters & Nuclear Star Clusters 

Figure [5] shows T e g (blue filled squares) and T c (red) for f g = 1 
and /dg.MW = 10 (open circles), 1 (filled circles), and 0.1 (open 
squares), versus inferred cluster stellar mass M, for clusters in M82 
(left panel; McCrady & Graham 2007) and for the nuclear star clus- 
ters of galaxies (right panel; Leigh et al. 2012). For the M82 clusters 
I computed E = M/AtvR 2 , using the values of the total mass, ve- 
locity dispersion, and crossing time quoted in McCrady & Graham 
(2007), while for the nuclear clusters I used M and R from Leigh 
et al. (2012). The bolometric flux was then assumed to be given by 
equation l[2}, assuming = 1500 Lq/M©. This procedure yields 
r e ff (blue squares) at cluster birth, assuming that the entire stellar 
mass formed in less than ~ 3 Myr. The optical depth was then es- 
timated using equation {8]l for the assumed /d g , mw, and equation 
lf5J was solved for T c , giving the red points. 

The separation in temperature between T c s and T c indicates 
that for f g — 1 and /d g ,MW > 0.1, most of the clusters shown 
in Figure[5]are optically-thick on the scale of their half-light radii. 
Since T c measures the minimum temperature in the optically-thick 
regions of protoplanetary disks embedded in these clusters during 
formation, we see that T c > Ti co in many systems, thus likely sup- 
pressing giant planet formation by core accretion or gravitational 
instability. 



5.2 Local Clusters & Globular Clusters 

In Figure [6] I plot temperature versus current stellar mass of both 
young and globular clusters of the Milky Way, Fornax, LMC, and 
SMC (McLaughlin & van der Marel 2005). The black open cir- 
cles show the current V^-band effective temperature computed us- 
ing T ff, nc ,w = (Lv /47ri?^crsB) 1,/4 , where Rh is the half-light 
radius and Lv is the V-band luminosity of the cluster within Rh. 
The blue points show the zero-age main sequence (ZAMS) effec- 
tive temperature computed using T cff = (*E/ctsb) 1/4 , where 
E = M/iirRh, M — TyLv is the total stellar mass within Rh, 
Tv is the V-band mass-to-light ratio, and ^ = 1500 L0/M© is 
the assumed ZAMS light-to-mass ratio. The blue points approxi- 
mate the effective temperature for the cluster in its first ~ 3 Myr 
and are directly analogous to the blue squares in Figure [5]for M82 
and the nuclear star clusters. On timescales longer than ~ 3 Myr, 
T e ff oc t^ 1 / 3 , as shown in Figure[Tj 



In order to estimate which clusters may have been optically- 
thick to the re-radiated FIR emission from dust, I assumed that 
f g = 1 and that the dust-to-gas ratio relative to the Galactic value 
is given by /d g , mw = 10' Fe ^ H ' (see eqs.[6]&[7]l. Using this pre- 
scription for /d g . mw, I find that only six clusters in the sample are 
optically thick on the scale of their half-light radii (all are globu- 
lars, labeled in red), and that they are only marginally so. By the 
same calculation used in the panels of Figure [5] for T c , only NGC 
6440, 6338, and 6441 would have T c ~ Ti cc , and although 47 Tuc, 
NGC 6266 (M62), and 6526 are optically thick at R h ,T c < Ti cc 
(marginally) by that calculation on the scale of Rh. 

Motivated by the excellent radial profiles presented in 
McLaughlin & van der Marel (2005), and the fairly strong changes 
in E(r) and p(r) within each cluster, I was motivated to go a step 
further and attempt to calculate the fraction of the total mass that 
might reasonably have T c > Ti cc . In order to solve this prob- 
lem fully, one would need to solve the radiation transport prob- 
lem in spherical symmetry, with distributed sources, of emissivity 
j(r) — <Pp(r), subject to the free function of the gas/dust density 
profile, p g (r). Moreover, one would need to infer the initial stellar 
density profile at ZAMS from the density profile we observe to- 
day. Because these steps are highly uncertain, I made the following 
crude approximation. I assume that at every radial position r, the 
total flux carried by that shell is the volume integral of the sources 
< r: 

f-r 

F(r) = r~ 2 j(r')r' 2 dr'. (30) 
Jo 

At that same position, the total overlying optical depth is 

roc 

r(r) = / K R f gP {r')dr'. (31) 

J r 

For the required density and emissivity profiles, I used the power- 
law (7-) models presented by McLaughlin & van der Marel (2005) 
for each cluster, which yield the three-dimensional emissivity pro- 
file 

j(r) = *p [1 + (r/ro) 2 ] ~ 7/2 , (32) 

where po and 7 are provided for each cluster by the authors. I then 
assumed f g = l and / dg , mw = 10' Fc//H ' , and solved equation l[5]l 
for T c (r) using F(r) = asBlcff and r(r) in equations J30b and 
< |3 lb - This procedure accounts for the fact that both the total flux 
and optical depth vary strongly with r, and allows one to make an 
estimate of the central temperature of the cluster where F(r) may 
be very different from F(Rh), and r(r) may be much larger than 
r(R h ). 

The solid red dots in Figure|6]show the resulting central clus- 
ter temperature. The small red number shows the fraction of the 
mass that has temperature larger than Ti cc ~ 150 K, M(T > 
T\ ce )/M. This calculation thus implies that only ~ 0.1 — 0.4 of the 
total mass of these clusters exceeds the ice-line temperature. Gener- 
ically, it is the central regions, where the flux is still very high, and 
the overlying optical depth is very large, where T c (r) > T\ ce ob- 
tains. I emphasize that this procedure is very crude since it (1) is 
subject to the uncertainties of a free function, p g (r), which I have 
specified to be f g p(r), (2) does not in fact solve the transport equa- 
tion with the non-gray opacity of dust, and (3) assumes that p(r) is 
unchanged from its initial distribution. Nevertheless, this procedure 
allows one to explore the obvious change in F(r) and r(r), which 
undoubtedly lead to temperature profiles in clusters at formation. 

I discuss a number of specific systems further in Section|6] but 
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Figure 5. Effective temperature at fonnation T c g (blue rilled squares) and T c , assuming 'if = 1500Lq/Mq and fa* MWi fg = ( 10, 1) (open red circles), 
(1, 1) (filled red circles), and (0.1, 1) (open red squares), for M82 super star clusters (left; McCrady & Graham 2007) and the central star clusters of galaxies 
(Leigh et al. 2012). Protoplanetary disks embedded in optically-thin clusters reach T c fj (blue). Disks in optically-thick clusters (red) reach T c during the first 
~3 - 4Myr. 



here note that none of the clusters highlighted in red, with M(T > 
Ti ce )/M > are thought to have undergone core collapse since 
birth. 



5.3 Galaxies & AGN 

The formalism discussed so far for star clusters can be general- 
ized to continuously star-forming galaxies by taking F — £E*c 2 , 
where E* is the star formation rate per unit area, and e is an IMF- 
dependent constant of order 5 x 10~ 4 — lCP' 5 for standard IMFs. 
The effective and central temperatures of galaxies in formation 
are then (compare with eqs. [2]and[9] see the related discussion in 
Thompson, Quataert, & Murray 2005) 

T cff ~102k(£_ 3 E*, 3 ) 1/4 (33) 

and 

T ~ 170 K [/dg, MW £-3 E g , 4 E*, 3 ] , (34) 

respectively, where s _ 3 = e/10" 3 , E*, 3 = E*/10 3 M Q yr" 1 
kpc -2 , and E g> 4 = Eg/10 4 Mq pc -2 . Equation d 34| > assumes 
t > 1 and equation l[6j[3 

As an application, in Figure [7] I show temperature as a func- 
tion of stellar mass estimated for the compact high-z galaxies from 
the sample of Van Dokkum et al. (2008) and Kriek et al. (2008). 
The filled black triangles show the current effective temperature 



10 Using e = 5 X 10~ 4 decreases eqs. g3] and [34] to 86K and WOK, 
respectively. 



Tcff, now = (M/T47ri? 2 ff o-sb) 1//4 , where T is the current mass- 
to-light ratio of the stellar population found by Kriek et al. (2008), 
and M and R e g are the total stellar mass and effective radius. I 
then compute the dynamical timescale tdyn ~ (Gp)~ 1//2 for each 
galaxy and assume that the star formation timescale is tsF = 10 x 
(open squares), 50 x (filled circles), 100 x tdyn (open circles). The 
flux at formation is then F = eEc 2 /fsF, where I approximate 
E = M/4tvRI s . The blue points show T cff = (F/a SB ) 1/4 , the 
effective temperature during formation, while the red points show 
the solution to equation ((5) assuming that E 9 = E (f g = 1) and 
/d g ,MW = 1, which approximates the estimate of equation ( 1341 at 
high optical depths. 

Although highly simplified, these estimates imply that all of 
these systems were > 30 K at formation and that many were 
> 70 K. And ~ 6 of the systems shown here likely attained tem- 
peratures in excess of Ti ce , almost certainly suppressing planet for- 
mation by core accretion or gravitational instability. 

Similar estimates can be made for actively star forming galax- 
ies. Equations l |33t and d34| l show that it is the most gas-rich and 
highly star forming systems that are likely to suppress planet for- 
mation. These include the inner regions of Arp 220 and other lo- 
cal ULIRGs, high-z submillimeter selected sources, and the highly 
star-forming clumps within high-z galaxies like that in Q2346-BX 
482 (Genzel et al. 2008; Murray et al. 2010). One can use the 
Schmidt law as derived by Kennicutt (1998) to connect E* and 
E s in actively star-forming galaxies. Equation ((34) then allows an 
estimate of the critical gas surface density above which T exceeds 

E g , i ce ~ 2 x 10 4 M pc" 2 Tl'l 150 K £-° 3 4 . (35) 
Such surface densities are in fact achieved in the inner regions of 
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Figure 6. Effective temperature now inferred from V-band light T^ w v (black open circles), effective temperature at zero-age main sequence T of j (blue filled 

squares), assuming ^ = 1500 Lq /Mq, and central cluster temperature T c (see Section 1531 assuming / 9 = 1 and /d g mw = 10l FC//H ' (red) for clusters 
in the Milky Way, Fornax, LMC, and SMC (McLaughlin & Van Der Marel 2005). The fraction of the total mass above the ice line temperature in the 6 
optically-thick globular clusters M(T > Ti cc ) /M is noted. 47 Tuc is highlighted (dot+circle). Protoplanetary disks embedded in optically-thin clusters reach 
T c ff (blue), whereas disks embedded in optically-thick clusters (red) reach T c during their first ~ 3 — 4Myr. 



Arp 220 (Dowries & Solomon 1998). Matsushita et al. (2009) show 
that this starburst is optically-thick at 435 fim, thus justifying the 
estimate that r > 1, and Downes & Eckart (2007) find a true dust 
temperature of 170 K for the inner disk. Herschel observations pre- 
sented by Rangwala et al. (2011) show that the dust has optical 
depth of ~ 5 at 100 /im, and multiple molecular line diagnostics 
indicate the presence and pervasiveness of hot dust and gas (see 
also Ottetal. 2011). 

Again, equation d33 1 sets the fundamental temperature scale 
for the protoplanetary disks embedded in highly star-forming 
galaxies and starbursts. Since star formation surface density and 
gas/dust surface density are highly correlated through the Schmidt 



law, regions with > 10 3 M yr 1 kpc 2 will have r > 1 and 
the temperature of the system will be higher than the effective tem- 
perature by a factor of r 1 / 4 . One sees from these simple estimates 
that, although extreme, many systems may reach this limit. 

Since the compact high-z galaxies presented in Figure |7]may 
be the cores of present-day ellipticals (Hopkins et al. 2009), our ex- 
pectation for the prevalence of planets in these metal-rich systems 
should be adjusted accordingly strongly downward. Although rare 
in the local universe, the conditions producing Ei cc in equation d351 > 
were more common in the high-z universe, and the subsequent high 
temperatures and suppression of planet formation via gravitational 
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instability or core accretion should affect the present day planetary 
mass density. 

Finally, AGN activity can also heat the medium above 
Ti ce . For a BH luminosity of L = L 4 6i0 46 ergs s~\ T eff > 

Ti C e for R < 50 pc -I/ 4 g ^ice, lsoK- Assuming an isothermal 
sphere for the stellar mass distribution in the central regions of 
galaxies, this corresponds to M(< R) ~ a 2 R/G ~ 5 x 
10 8 M (o-/200km/s) 2 (i?/50pc), where a is the velocity dis- 
persion, about 1% of the total bulge mass. Supermassive BHs may 
undergo significant growth during rapid star formation episodes on 
multi-kpc scales. The planet population may thus be affected by 
both the episode of BH growth and the star formation itself. 



6 DISCUSSION 

6.1 Caveats, Complications, & Uncertainties 

Radiation Transport — Equation (O is highly approximate. It ap- 
plies only in the limit of a an overlying column of matter without 
sources, makes the Eddington approximation, and assumes both an 
infinite plane-parallel atmosphere and gray opacity. A real embed- 
ded star cluster has a distribution of sources (stars), whose light is 



reprocessed locally into the infrared by absorption, and which con- 
tributes to the local flux in the atmosphere. The radial distribution 
of sources j(r) = 9p(r), the distribution of dusty gas with respect 
to those sources, f g (r) = E 9 (r)/E(r), the variation in the IMF 
as a function of radius ^z(r), the variation in the dust-to-gas ratio 
/d g , mw(0 and thus the opacity, will all affect the radial run of 
temperature with optical depth, T(t). 

A detailed exploration of these effects is beyond the scope of 
this work, but I note here that the formalism of Hubeny (1990) for 
distributed sources in plane-parallel atmospheres can be general- 
ized to spherical systems with non-gray opacities (see, e.g., Whe- 
lan et al. 201 1). However, more detailed radiation transport models 
should be coupled to a self-consistent calculation of the dynamics 
of star-forming and disrupting clusters since f g (r) = E 9 (r)/E(r) 
would otherwise be a free function. In addition, since star forma- 
tion may be spatially correlated within individual star clusters, the 
radiation transport problem is likely to be further complicated by 
multi-dimensional effects (e.g., Krumholz et al. 2012), and, further, 
by the dynamics of merging sub-clumps. 

The Optical Depth & Turbulence — Recent calculations sug- 
gest that massive clusters may be disrupted by radiation pressure 
on dust grains, and this feedback mechanism has been suggested 
as a way to regulate star formation on galaxy scales (Murray et 
al. 2005; Thompson et al. 2005; Krumholz & Matzner 2010; Mur- 
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ray et al. 2010; Andrews & Thompson 2011; Hopkins et al. 2012). 
The parameter space over which this feedback mechanism may 
dominate largely overlaps with the parameter space discussed in 
this work (e.g., E > 10 3 M pc" 2 or E* > 10 3 M Q yr _1 kpc~ 2 ). 
If the gaseous medium is highly turbulent, and in particular, if the 
radiation field itself drives that turbulence, then the effective op- 
tical depth for the reprocessed IR radiation may be less than the 
simple estimate of equation {8} (see Murray et al. 2005; Hopkins 
et al. 2012). Recent simulations by Kramholz & Thompson (2012) 
of sustained radiation pressure driven Rayleigh-Taylor convection 
(in 2D) in fact suggest that in these environments the value of the 
effective optical depth is decreased by a factor of ~ 2 — 5, as low 
optical depth channels through the overlying medium are formed 
and the radiation field more efficiently escapes. In principle, this 
should decrease the temperature below the approximation of equa- 
tion (f5}, but this effect has yet to be quantified in the super star 
cluster environment. 

IMF — Observations imply that the light-to-mass ratio of a 
ZAMS stellar population may vary. The value of ^ I scale to 
throughout this work, I5OOL0/M0, is approximately appropriate 
for a Salpeter IMF from 1 — 100 M Q . More precisely, <3/(t = 
0) ~ 1400 and #(t = 3Myr) ~ 19OOL /M Q (see Fig. Hi. 
For limits of 0.1 — 100 Mq, unrealistic for the low-mass stellar 
field population (e.g., Kroupa et al. 1993; Kroupa 2001; Chabrier 
2003), <&{t = 3Myr) decreases to ~ 750. For a Kroupa IMF 
i£(t = 3Myr) ~ 12OOL Q /M . A somewhat more bottom-heavy 
IMF below 1 Mq, but not as steep as Salpeter, as observed in the 
Galactic bulge by Zoccali et al. (2000), would yield similar values 
for ^. Stolte et al. (2006) measure the present day mass function 
in the young compact Galactic star cluster NGC 3603 and find a 
slope of N oc M~°' 9±0 ' 15 from ~ 1 - 20 M . In Arches, Stole 
et al. (2005) find a similar slope at high masses, but with a signif- 
icant turnover at ~ 6 — 7M . Such a flatter than Salpeter slope 
at high masses, if extended to ~ 100 Mq, and taking Kroupa from 
0.1 — 1 Mq, would increase <]/ over our nominal value by a factor 
of ~ 1.5 - 2. 

The value of and its time dependence also change if the 
stellar population is rotating (Ekstrom et al. 2012). Levesque et 
al. (2012) (e.g., their Fig. 6) show that ^/ reaches a maximum 
at ~ 3.5 Myr, effectively extending the life of the massive stars 
that dominate the bolometric output of the stellar population by 
~ 0.5 — 1 Myr over the non-rotating population. At an age of 5 Myr, 
the rotating population is ~ 1.6 times brighter for the same stellar 
mass, and at 10 Myr it is ~ 1.3 times brighter. Rotation would thus 
increase the temperatures of the protoplanetary disks of star clus- 
ters, all else being equal. 

The normalization, shape, and high-mass slope of the IMF 
may also change as a function of cluster mass and surface den- 
sity. Murray (2010) presents an argument that it is precisely the 
FIR optically-thick star clusters with the highest values of E and 
E 9 that will deviate most strongly from a standard IMF, with larger 
ty. Such a correlation — 9 oc E or M to some power — would 
strongly increase T e g and T c with respect to equations $2$ and ® in 
the highest surface density systems. As an example, super star clus- 
ter M82-F was argued by Smith & Gallagher (2001) to be strongly 
inconsistent with a standard Kroupa IMF, and requires a low-mass 
cutoff of ~ 2 — 3M . Such an IMF would increase for the 
ZAMS population, increasing the inferred temperature in the first 
few Myr of evolution. A strongly bottom-heavy IMF, as has been 
claimed in elliptical galaxies (Van Dokkum & Conroy 2010, 201 1; 



Wegner et al. 2012), and which may be a function of velocity dis- 
persion (Spiniello et al. 2012), would go the other way, decreasing 
by a factor of ~ 2 or more. 

Finally, there is evidence for a bottom-heavy IMF in ellip- 
tical galaxies (Van Dokkum & Conroy 2010, 2011; Wegner et 
al. 2012), which may be a function of velocity dispersion (Spiniello 
et al. 2012), and which would decrease the temperature by decreas- 
ing ^ relative to the fiducial value used throughout this work. 

Time Dependence — The central temperature of the star clus- 
ter is likely to change as a function of time, as the gas is used up, 
the stellar population grows, the stellar population evolves dynami- 
cally, and as the cluster is disrupted, ejecting the remaining gas and 
perhaps unbinding or expanding the stellar distribution. A more re- 
alistic treatment would couple the evolution of T throughout the 
cluster with a dynamical calculation of cluster disruption, as in the 
simplified models of Murray et al. (2010), or in the (much lower 
mass) cases of numerical star cluster formation with realistic ther- 
modynamics (e.g., Bate 2010, 2011, 2012; Krumholz et al. 2012). 

In any case, since the timescale for the embedded cluster phase 
is likely to be a multiple of the dynamical timescale, the average 
inner cluster temperature may decrease rapidly from T c to T e g 
as the dusty column trapping the radiation disperses. Since this 
timescale is of order the timescale required for planet formation 
~ Myr, planet formation may be fully quenched by the fact that the 
medium is above Ti cc , or merely delayed for ~Myr until cooler 
climes prevail, T c < Ti cc and giant planets form. 

A related point is that the planet formation timescale is short- 
ened in the cluster environment because of the increased heating of 
the surface layers by the cluster irradiation (Section [4.3 1 , radioac- 
tivity of the (potentially) large-scale starburst environment (Lacki 
2012a), the very large cosmic ray ionization rates and gamma- 
ray fluxes expected in those environments (Lacki et al. 2010; Pa- 
padopolous 2010; Papadopolous et al. 2011; Lacki 2012b), or the 
UV irradiation from nearby massive stars (see Section[T). All act to 
increase the size of the active accreting zone in layered disks (Gam- 
mie 2001), and hence speed both the accretion of the disk onto the 
host star, and potentially the migration of planets in formation. 



6.2 RV & Transit Surveys in Galactic Globular Clusters & 
Open Clusters 

As discussed in Section loTTI there are substantial uncertainties in 
the estimates of T e g (blue), and especially T c (red) in Figures [5] [6] 
and|7] Even so, they make it clear that globular clusters, star clus- 
ters in the local universe, compact high- 2 galaxies, and nuclear star 
clusters were born very hot, and that many of these systems were 
likely optically thick at birth, leading to T c above the ice line tem- 
perature Ti ce . Given the discussion of Section [4] one would thus 
expect that these systems should not have been able to form gas- 
and ice-giants by core-accretion during their first ~ 1 — 10 Myr of 
evolution, that formation of similar systems by gravitational insta- 
bility would have been prohibited in all but the most massive disks, 
and that there should not have been a delivery mechanism for water 
to terrestrial planets on these short timescales. 

In Figure I highlight 47 Tucanae ([Fe/H] = -0.76; 
McLaughlin & van der Marel 2005), which has been the subject of 
an intensive observational campaign to find giant planets via tran- 
sits (Gilliland et al. 2000; Weldrake et al. 2005). Using our nominal 
numbers for 47 Tuc — / dg , M w = 10 [Fc/H] and f g = 1 — I find 
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that ~ 0.1 of the mass had T c > Ti cc . Changes of a factor of 2 in 
f g , ty, or the effective radius for the stars change this fraction from 
~ 0.1 to~ 0- 1. 

Given these substantial uncertainties, it is worth considering 
the implications of the hypothesis that the lack of planets observed 
in 47 Tuc is a result of a high temperature at birth. There are several 
possibilities, depending on whether or not planets form predomi- 
nantly by core accretion or gravitational instability, and what the 
true frequency of giant planets is as a function of metallicity, in- 
dependent of birth environment. If planets form predominantly by 
core accretion, and if the lack of giant planets in 47 Tuc in the sur- 
vey of Gilliland et al. (2000) is indeed entirely because T c > Ti CB in 
its first Myrs of evolution, then I predict that globular clusters with 
similar metallicity, but lower E will harbor giant planets if they 
had T c < Tice. For this reason, I highlight NGC 6366 in Figure 
which has a metallicity comparable to 47 Tuc ([Fe/H] = —0.82; 
McLaughlin & van der Marel 2005), but an average stellar surface 
density that is ~ 10 times smaller. NGC 6366 is the only relatively 
high-metallicity globular cluster with E substantially below 47 Tuc 
cataloged in McLaughlin & van der Marel (2005), and it is thus 
an ideal case for testing planet formation as a function of surface 
density at nearly fixed metallicity. 

A similarly important test is the globular cluster NGC 6121 
(M4), where a 2.5 Mj up i t cr mass planet has been inferred by Sig- 
urdsson et al. (2003). NGC 6121 has [Fe/H] = -1.2, and a simple 
extrapolation of Fischer & Valenti's correlation between planet fre- 
quency and metallicity would suggest that it should harbor very few 
giant planets. Is this planet's existence in conflict with our hypothe- 
sis for 47 Tuc? No. As shown in Figure|6] the effective temperature 
of NGC 6121 was likely substantially below Ti ce at birth as a result 
of its low metallicity and surface density with respect to 47 Tuc. 

The planet in NGC 6121 can also be thought of as a test for 
gravitational instability, which would then have to operate in a ra- 
diation environment of T e g « 40 K, implying a massive protoplan- 
etary disk (see Section l4~2l . 

Several additional systems are worth discussing in this con- 
text. In particular, the more metal-rich globular clusters of the 
Galaxy listed in McLaughlin & van der Marel (2005) are NGC 
6440, 6441, and 6388 with [Fe/H] = -0.34, -0.53, -0.60, re- 
spectively. These are the hottest three systems T c significantly 
above Ti ce , shown in red in Figure [6] Therefore, these systems 
also present a test of the core-accretion scenario, and this paper. 
If 47 Tuc indeed formed no (or few) gas- and ice-giants because 
T c > Tice, then NGC 6440, 6441, and 6388 should be similarly 
devoid. 

I emphasize that such a set of tests would be complicated by 
the many effects that may inhibit planet formation in dense clusters 
(see Section [TJ. Nevertheless, comparing the planet population at 
different metallicity and at different stellar surface density is one 
of the few ways of understanding these joint dependencies. 

Finally, note that one of the most compact and metal-rich of 
the Galactic open clusters is NGC 6791, which has been the subject 
of a number of transit surveys (Section [T}, has T c ff(ZAMS) ~ 
20 K and r R ~ 0.04, using M ~ 5000 M Q , R ~ 5 pc, / dg , mw = 
10 [Fc/H] , [Fe/H] = 0.3, and taking f g = 1 (Platais et al. 2011 
and references therein), and thus its planet population should not 
be strongly affected by the effects described in this work. 



6.3 Observations of Protoplanetary Disks in Dense Stellar 
Environments 

The simple calculations shown in Figure [3] and the discussion of 
Section[3]indicate that the vertical structure of protoplanetary disks 
in dense stellar environments should be strongly modified from that 
derived for passive or active disks without external irradiation. The 
work of Stolte et al. (2004), (2006), and (2010) and others use NIR 
excess to investigate the properties and incidence of disks within 
clusters like NGC 3603 and Arches. More detailed calculations of 
disks with self-consistent thermodynamics and structure are war- 
ranted for observations now and coming (e.g., JWST) of stars and 
their disks in the most massive and dense Galactic clusters. 

As an example, using the data collected in Table 2 from 
Portegies-Zwart et al. (2010), Arches has T eff ~ 100 K and T c ~ 
180 K for f g = /d g , mw = I- NGC 3603 has lower temperature, 
and is not optically-thick for f g = /d g , mw = 1, but still has 
T c ff ~ 70 K. Trumpler 14 and Westerland 1 have similar stellar 
surface densities and thus similar temperatures to NGC 3603. Fig- 
ure [3] and the scalings of Section [3] imply that the protoplanetary 
disks in these systems should be strongly modified by their radia- 
tion environment. 



6.4 The Maximum Temperature of Dense Stellar Systems 

The maximum stellar surface density observed in star clusters, el- 
liptical galaxies, globular clusters, and other dense stellar systems 
is (Hopkins et al. 2010) 

E max ~ 10 5 M Q pc' 2 . (36) 

If this stellar surface density was formed in a single burst, F max = 

1.5 x 10 14 #isoo L Q kpc" 2 , T e 'g ax ~ 180 K #^ , and r^ lax ~ 
100 /dg. mw fg- Taking the maximum dust-to-gas ratio observed 
in galaxies to be /d g , mw ~ 10 (Munoz-Mateos et al. 2009) im- 
plies that the maximum temperature it might be possible to achieve 
within a star cluster is 

T max ~ 10 3 K (/dg.Mw f a #i5oo/10) 1/4 . (37) 

Compare with Figure [2] This temperature is equivalent to the ef- 
fective insolation at a distance of ~ 0.5 AU around star of Solar 
luminosity, and it is perhaps the largest value of the IR cluster in- 
solation that can be obtained in the central regions of dense stellar 
systems. Although T max exceeds the condensation temperatures of 
magnetite (~ 370 K), troilite (~ 700 K), and approaches the con- 
densations temperatures of feldspar and forsterite this value is still 
substantially below the sublimation temperatures of the primary re- 
fractory condensates, and so it should not be possible to directly 
inhibit the formation of rocky planets. 



6.5 The Fraction of All Star Formation Occurring above a 
Given Stellar Surface Density 

If massive dense star clusters and galaxies in formation inhibit the 
formation of gas and ice giant planets by core-accretion, as implied 
by Sections|4]and[5] then it is worth asking about the fraction of all 
star formation that has occurred above the critical surface density 
Eicc, such that T c > Ti ce , as estimated in equation d21t . Such a 
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calculation requires several ingredients. Equation (|9) implies that 

dlnT c = idlnE+ idlnZ+ idln*, (38) 

where I have assumed that f 9 ~ const, and Z is the metallicity nor- 
malized such that /d g , mw oc Z, which enters into the dust opacity 
per unit gram of gas in equation ((6}. Over the history of most of the 
star formation in the universe, one expects Z to vary by ~ 2 — 3 dex. 
The total range in <]/ is likely smaller, only varying by (at most) a 
factor of ~ 3 — 5, and it may be approximately constant. The total 
range in E is quite large, spanning from ~ 10 2 — 10 5 Mq pc~ 2 . 
Naively, then, the fraction of all stars formed above Ei cc depends 
primarily on the cluster mass function dN/dM, but also dN/dR 
and dN/dZ, with a much smaller dependence on dN/dty. 

Typical functions for dN/dM for clusters yield dN/dM oc 

M~ 2 , but with a cutoff at high masses that is dependent on the 

galaxy observed. It is natural to take the maximum mass to be a 

fraction of the Toomre mass in a marginally stable galactic disk 

of gas surface density E 9 and scale-height h, 7r/i 2 E 9 , but the star 

formation efficiency as a function of the collapsed gas cloud may 

vary as a function of its mass. Similarly, the relationship between 

cluster mass or richness and its size is likely not a simple power-law 

over the range from M ~ 10 2 - 10 7 M Q (e.g., see discussions in 

Murray 2009; Murray et al. 2010; Adams 2010). In particular, R oc 
M o.3-o.4 for M < 1Q 4 Mq (Lada & Lada 2003) and R a M o for 

10 4 < M < 10 6 M with R ~ 2 pc. For larger mass clusters the 
relationship apparently steepens, with R oc M ' 6 (M > 10 6 Mq) 
(Walcher et al. 2005; Murray 2009). 

Given these various uncertainties and unknowns I take a sim- 
ple tack and make an estimate for an average galaxy that domi- 
nates the star formation rate of the universe at redshift ~ 1 — 2. 
Such star-forming galaxies have gas surface densities ~ 10 — 100 
times the Galaxy, with commensurately higher star formation rates 
per unit area. Their metallicities are likely sub-Solar on average 
by a factor of a few. If one assumes that a typical cluster size 
is that dN/dM oc M~ 2 , then the fraction of all star formation 
occuring with T c > Ti cc ~ 150 K is just f(T c > Ti ce ) = 
ln(M max /M crit )/ ln(M max /M min ) where 

M crit -2.5 x 10 s M ( J R/2pc) 2 (/ dg ,Mw/ 9 *i5Oo)" 1/2 (39) 

is the critical cluster mass above which T c > Ti ce (see eq.!21b. and 
M m i n and M max are the minimum and maximum cluster masses, 
respectively. Taking R = 1,2, and 4 pc, /d g , mw = fg = *isoo = 
1, M mln = 10 M Q , and M max = 1O 7 M , / - 37%, 27%, and 
17% of stars are born in clusters with T c > Ti ce , respectively. For 
/dg, mw = 0.1, these fractions decrease to ~ 28%, 18%, and 8%, 
respectively. 

In this context, it is worth noting the recent work of Gieles 
et al. (2010) who provide an explanation for the Al — R rela- 
tion for stellar systems above — 10 4 Mq (see their Fig. 3). They 
show that dynamical processes increase the half-light radii (R) of 
clusters of Gyr timescales. Their initial configuration starts with 
M ~ 10 6 Mq (i?/pc) 2 , a constant surface density at all masses 
equivalent to E ~ 8 x 10 4 Mq pc -2 and an effective temperature at 
zero-age main sequence of T e g ~ 170 K from equation J2j, Since 
the timescale for significant dynamical evolution in these simula- 
tions is lOMyr or longer I expect systems that start with this value 
for E to have T c > Ti ce throughout the crucial planet formation 
epoch. If such an initial function for the M — R relation was in 
fact generically valid for M > 10 4 Mq, the arguments presented 
in this paper would imply that more than 50% of all stars were bom 



in an environment hostile to formation of gas and ice giant planets 
via core accretion or gravitational instability. 

Finally, in a study of the structure of elliptical galaxies, and 
compact galaxies at high redshift, Hopkins et al. (2009) estimates 
that the total mass in present-day ellipticals with stellar surface den- 
sities in excess of Ei ce is — 10 — 20%. 



7 CONCLUSIONS 

The basic conclusions of this work are as follows: 

1. Massive star clusters with standard IMFs that form in a sin- 
gle burst of star formation are hot for a time comparable to the 
planet formation timescale, reaching fluxes equivalent to an effec- 
tive temperature of — 50 — 150 K for — 3 — 4 Myr (FigureQ] eq.|2j- 

2. Embedded dusty star clusters reach a temperature larger 
than T e g because they are optically-thick to their own reradiated 
FIR emission. For surface densities above ~ 5 x 10 3 Mq pc~ 2 , 
Galactic dust-to-gas ratio, and gas fractions at birth of order unity, 
the characteristic temperature of star clusters exceeds the ice line 
temperature (Figure|2] eqs.l9]|2"TV See the discussion of uncertain- 
ties and caveats in Section [6~T| 

3. Protoplanetary disks in such clusters are strongly isotropi- 
cally irradiated, and this cluster irradiation dominates the thermo- 
dynamics of both passive and active disks around Solar-type stars 
at semi-major axes beyond — 1 — 5 AU, changing the disk structure 
(Figured eqs. U2ll4l and[T8l&[T9l. 

4. In the most massive dense star clusters, the temperature of 
protoplanetary disks should exceed the temperature for water ice 
condensation at — 150 — 170 K (Figures[5]&[6]l while the cluster is 
embedded. This decreases the total amount of condensable material 
by a factor of ~ 2 — 5, prohibiting formation of gas- and ice-giant 
planets if they form by core accretion (Section r4.il . and prohibiting 
habitability. 

5. Planet formation by gravitational instability will also 
be suppressed in protoplanetary disks within clusters because 
Toomre's Q > 1. Typical disk masses that could become gravi- 
tationally unstable are required to exceed ~ 50 — 100 times the 
Minimum Mass Solar Nebula (Figuregl Section|4~2l eqs. 122124b . 

6. In these highly-irradiated environments, the overall accre- 
tion rate should be enhanced, leading to shorter disk lifetimes and 
less time to build planets. The migration rates of planets in gas or 
planetesimal disks are also likely to be influenced (Section [4.3b . 

7. Even accounting for their low metallicities, 47 Tuc and sev- 
eral other Galactic globular clusters plausibly exceeded the limiting 
surface density for their central regions to exceed Ti cc , and this may 
be a factor in explaining why 47 Tuc appears to have a paucity of 
hot Jupiters (Figure [6] Section 15.2b . A test of this scenario would 
be an RV/transit survey of NGC 6366; with a metallicity similar to 
47 Tuc, but a much lower surface density, it should have abundant 
giant planets. Conversely, the relatively more metal-rich globulars 
NGC 6440, 6441, or 6388, which had even higher temperatures 
than 47 Tuc at birth should be devoid of planets (Figure [6] Section 
16.21 related discussion of NGC 6121). A mutual comparison be- 
tween planet populations in clusters of different metallicities, but 
at fixed surface density, and of different surface densities, but at 
fixed metallicity, would be very valuable. 
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8. Rapidly star-forming galaxies also maintain very high tem- 
peratures. Those with average star formation rate surface densities 
> 10 3 Mq kpc -2 yr _1 have central temperatures that approach 
Tico, and should thus globally suppress planet formation (eq. 1341. 
AGN activity might too. Simple estimates suggest that the quies- 
cent compact high-z galaxies plausibly attained temperatures ex- 
ceeding Tico during formation (Figure [7}. The idea that these sys- 
tems might be completely devoid of giant planets and habitable 
planets is startling. Because these systems may be the cores of to- 
day's elliptical galaxies, and because a component of the Galac- 
tic bulge and M31's bugle may have formed in such an extreme 
starburst, one generically expects future surveys of the planet pop- 
ulations in these environments (e.g., via microlensing) to show a 
paucity of giant planets per star, even though (and, indeed, because) 
these systems are metal-rich. 

9. Finally, although the estimate is uncertain, between ~ 5 — 
50% of the z ~ stellar mass density formed in a star cluster 
with temperature exceeding Xi ce , and thus our expectations for the 
inventory of planets, and their incidence per star as a function of 
birth environment and metallicity may be significantly impacted. 
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